function C = diffEq(t,a,b,CF,C0)

tmp = (1 - b*t)./sqrt(4*a*t);
tmp(1) = 1e10; 
C = (CF - C0)*(1 - erf(tmp))/2 + C0;
